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Abstract 

CN ; 

t— ( \ The ground state of the four-chain Heisenberg ladder model is numerically investigated. Hidden- 
^ order correlations suitable for the system are introduced and calculated with an emphasis on the 
q\ • spatially isotropic point, where a corresponding material exists. The existence of a long-range 
O i hidden correlation indicates formation of a short-range RVB state in the case of the antiferro- 
q . magnetic inter-chain coupling. A transition between the phase of the ferromagnetic inter-chain 
| coupling and that of the antiferromagnetic one is discussed. 

c3 
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i ■ 1 Introduction 

O In recent years, considerable attention has been devoted to ground-state properties of low- dimensional 
H quantum systems. The interest was greatly stimulated by the discovery of the high-temperature 
[ superconductivity, and by subsequent studies which revealed the importance of the Cu0 2 plane 
^ ■ of high-T c materials. However, complete understanding of the high-T c mechanism still appears 
■ to be beyond our reach. A possible new approach to the high-T c mechanism was provided by a 
recent experiment on the novel material Sr2 nicg _2Cu2n lcg 04 nicg _2. This material has periodic line 
defects in the Cu02 plane, and thus is composed of ladders with rii eg legs interacting weakly with 
each other. The limit rii eg — > oo yields the high-T c superconductors. 

At half-filling this material may be described well by the antiferromagnetic Heisenberg model 
defined on a ladder: 

"leg L "leg— 1 L 

t~c = j y, y, Sn ■ Sn + i + j' y y sm ■ Si+ij. (i) 

i=i i=i i=i i=i 

(Here denotes the 5 = 1/2 spin at the zth site of the Zth chain. We put J = 1 hereafter.) 
It was reported, indeed, that the double chain (ni eg = 2) with J ~ J' develops the paring 
correlation upon doping. This report was followed by the suggestion ||, |5|, || that a set of ladders 
interacting weakly with frustration might show the superconductivity. Thus, characterization 
of the ground state of the ladder models can be an essential step to understanding the high- 
temperature superconductivity. 

It has been conjectured0, [5], |J for J = J' that the ground state of the model (|l|) with even rii cg 
is spin liquid with the energy gap, while that with odd n\ eg is critical, or gapless. This conjecture is 
supported by a theorem, numerical calculations up to ri\ eg = 4,|| |], |] and a scaling theory. [|10|1 



In order to explain this remarkable conjecture intuitively, White et a/.0 proposed a resonating- 
valence-bond (RVB) picture of the ground state of the antiferromagnetic ladder models. This was 
followed by a proposal [|Tl| , [12 , [T~3|1 of a hidden-order correlation for the RVB ground state of the 
double chain: 



O 



RVB \ 



(\i = ({S" u + Sl t+l )e i ^^ +s ^\St 1 + S 



2,j+l- 



The existence of the long-range RVB correlation supported [|TT], |Tj| the RVB picture for ri\ eg = 2. 

The present RVB correlation (0) was also shown ||Tlf to be very useful in discussing another 
issue, namely the criticality of the model ([!]). Suppose that we change the value of the coupling 
across the ladder, J' (the inter-chain coupling). Then the above conjecture arises the following 
question. The model is known to be critical, or gap less [14] for J' = 0, whereas the model has the 
gap for J' = J according to the above conjecture. How does the energy gap emerge as J' changes? 



In the previous paper |LT| we presented numerical calculations of the RVB correlations as well as 



another hidden-order correlation, namely the string correlation: |TS] 



0, 



string 



(3) 



We showed |TTJ quite clearly for n\ eg = 2 that the energy gap appears as AE ~ J' u with u = 1. 



This conclusion is consistent with other studies. 10, 16, 17, 18 



Hidden correlations were first introduced to investigate the ground state of the integer-,!} anti- 
ferromagnetic Heisenberg chains. The hidden correlations relevant to the model are the following 



string correlations, [19, 20, 21, 22, 23, 24, 25 



X, Z) 



(4) 



where {Si} denote the S = 2 spins. The correlations successfully detected the valence-bond- 
solid (VBS) structure p6| ^7], ^ of the ground state of the integer-spin chains as well as their 
criticality. 0, g2| 

Now that we know the usefulness of the RVB correlation in the case n\ eg = 2, it is a challenging 
problem to generalize the definition (^) to higher n\ eg , and to explore the ground state and the 
criticality of general n\ eg ladder models. Here we consider the ground state of the ladder model 
with four legs. The Hamiltonian is given by (|1|) with n\ eg = 4. We generalize the definitions (|2|) 
and (0) for the four-leg ladder, and show the RVB structure of the ground state in terms of the 
correlations. We also confirm the prediction of the scaling theory [|10| that the ground-state phase 
diagram consists of two disordered phases covering J' > and J' < 0, respectively, and the critical 
point between them, J' c = 0. This phase diagram is the same as in the two-leg case. (16], [L7|, 18 

The ground-state properties are known in some cases. In the limit J' —>■ oo, the model is 
decoupled to independent four-spin rungs. The ground state is given by the product of singlets 
formed on the rungs. f29| Therefore, the energy gap is given by AE(J') ~ J' in this strong-coupling 
region. At the isotropic point J' = 1, which is of experimental interest, @ the energy gap AE 
and the correlation length £ have been numerically estimated as AE = 0.190 and £ = 5 ~ 
6.0 At J' = the model is reduced to four S = | Heisenberg antiferromagnetic spin chains. 
Hence, it is critical as noted above. In the limit J' — > — oo, the model converges to the S = 2 
Heisenberg antiferromagnetic chain, which is massive according to the Haldane conjecture. [30, pT| 

2 has been confirmed numerically. [11, 25, 32, 33, [34|] 



The Haldane conjecture for S 

The present paper is organized as follows. We introduce hidden correlations in the next section, 
and present a naive discussion on the existence of the hidden correlations in terms of the RVB 
argument. In section 3, we show numerical results of the hidden correlations with an emphasis 
on the case J' = 1, where the corresponding material is available. We also present an argument 
on the existence of the RVB structure by considering two double-chain ladders coupled weakly. 



This argument intuitively shows the nature of the RVB ground state. In section 4, we elaborate 
on the criticality of the model in terms of the hidden correlations. In the last section, we give a 
summary of the present paper. 



2 RVB state and hidden correlations 

In this section, we introduce the hidden correlations for the present four-leg-ladder model ([!]). 
They are defined so that they can detect expected RVB patterns. 

We present schematic drawings of the expected RVB patterns in Fig. 1 (a) for J' > and 
(b) for J' < 0, respectively. Each RVB pattern is arranged as follows. First for J' > 0, White 
et ct/.0 proposed that the ground state is dominated by the RVB pattern with vertical singlets 
and horizontal singlets. The state is stabilized by the resonance between the configuration of two 
adjacent vertical singlets and the configuration of two horizontal singlets; see Fig. 1 (a). This 
RVB picture was reported^ to be very useful for understanding several features of the ladder 
models. Second, the expected RVB pattern for J' < is given as follows. In the limit J' — > — oo, 
the system ([I]) with n\ eg = 4 converges to the S = 2 Heisenberg antiferromagnetic chain. The 



ground state of the chain was described [[22|, ^3], ^4 well in terms of the S = 2 valence-bond-solid 



state. This VBS state is constructed in the following manner: Suppose that each spin with S = 2 
consists of four spins with S = |; Form a singlet using two S = | spin of neighbouring sites; 
Arrange the singlets so that each bond may have two singlets. All the S — ^ spins are thus 
connected each other with the short-range valence bonds consequently. We expect this state in 
the limit J' — > — oo. The resonance may arise among the different pairing patterns for finite and 
negative J'. Hence we have drawn the pattern which is shown in Fig. 1 (b). We expect that this 
pattern dominates in the whole region J' < 0. 

In order to detect the above RVB patterns, we define the following hidden correlations: 

05vb(0, \* = (i S ti + S W + S h + ^, +1 )e i ^rJ(^ fc +^ fc+1 +^ fc +^ fe+1 ) 

(S^j + 5*2 j +1 + 5*1 j + S^ j+l )J (5) 

and 

(S- +Sl i + Sl i + Si ,)). (6) 



We here generalized the hidden correlations (0) and (0), which have been applied to the double- 
chain ladder model. [ll|, 0, 0] The angle n appearing in (§) and (|]) is replaced by 9 in the present 



definitions @ and (§). 

We show in the present paper that the choice 9 = | is the most relevant to the four-chain 
ladder model. We also show that the correlation develops in the phase J' > while the 
corelation (|j) develops in the phase J' < 0. Let us explain the reasons of these facts briefly, before 
going into details in §3 and 4 

First, in the limit J' — > — oo, the model converges to the S = 2 Heisenberg antiferromagnetic 
chain, as is explained in the previous section. The correlation then is reduced to 

(SfeV^Hsj), (7) 

where {Sj} denote the S — 2 spins. This correlation with 9 — ^ is reported to remain finite in 
the limit \i — j\ — > oo,[23, [ZIJ as was mentioned in §1. The correlation (Bl), on the other hand, 



remains finite in the limit J' — > oo; see the next section. 



Second, as was done for the double-chain ladder model, 01 we can show the existence of hidden 



orders schematically as follows. Consider any spin configurations of the four-leg ladder, satisfying 
the RVB pattern of Fig. 1 (a). We present an example in Fig. 2 (a). We notice that the set of 
the numbers {Sf = Sf i + S^i + S^ { + S%j} in Fig. 2 (a) satisfies the condition 



£3 

k=i 



< 2. 



(8) 



for v.- v 



j. This is due to the formation of two singlets over two neighboring composite spins {S^}. 



It is easy to see that the quantity 



(9) 

remains finite for any % and j, when it is averaged over all configurations that satisfy the condition 
(|8D.||52"|, [23|j This explains that the RVB correlation (EH remains finite if the RVB configurations as 
Fig. 1 (a) are dominant in the ground state. The above argument applies to the RVB patterns 
as Fig. 1 (b) if we define the composite spin as {S- = + S£ i+ i + + S% i+1 }; see Fig. 2 
(b). We can expect that local destruction of the RVB patterns does not result in vanishing of the 



long-range correlation. 11 . 



3 RVB formation and the presence of the hidden order 

for J' > 

In this section,we investigate the ground state of the four-leg ladder by means of the hidden 
correlations (|5|) and (||) numerically. The RVB picture explained in the above section is confirmed. 
We show that the choice 9 — ~ is the most relevant to the present system. 



3.1 Hidden correlations in the limit J' 



oo 



In this subsection, we concentrate on the four- leg ladder in the limit J' —>■ oo. The limiting point 
is expected to be the fixed point of the phase J' > O.fLOf. Hence, characteristics at the point may 
be relevant to those of the whole region J' > 0. 

In this limit J' — > oo, the Hamiltonian is reduced to a set of independent rungs. The ground 
state is given by the direct product of singlets that are formed on the rungs. We can calculate the 
long-range limit of the hidden correlations (^) and (|6|) exactly. Though the ground state is spin 
liquid. the hidden correlation (]5|) is long-ranged. The correlation Of-yB^ °°) is plotted 

against 9 in Fig. 3. The correlation has the maximum around 9 = n/2, while the correlations at 
9 = and tt are both vanishing. These behaviors are intrinsic to the whole region of the phase 
J' > as we see below. It is apparent, on the other hand, that the correlation M) is of short 



range, O. 



string \ 



OO 



0. 



3.2 The hidden orders at J' = 1 

In this subsection, we concentrate on the system at J' = 1, where the interaction is spatially 
isotropic. This isotropic system is of experimental interest. @. A strong-inter-chain-coupling 
expansion starts to fail at this point. 



We employed the density-matrix renormalization-group 
method |35|, |36| in order to treat large systems approximately. We conclude that ground state 
properties are indeed consistent with the RVB picture given in the previous section and the recent 



proposal that the point J' — >• oo is the fixed point of the phase J' > O.flTO 



First of all, we show the precision of the present renormalization-group calculation. In Fig. 
4 we plotted the relative error of the ground state energy with J' = 1 and L = 6 against the 



approximate level m. The parameter m is the number of states kept; [[5^, [SB] in this method, we 
treat only m states in the course of the renormalization. It should be noted that the precision 
of the correlations is worse than that of the ground-state energy. The reason may be as follows. 
Because the Hilbert space is restricted to the m states, this renormalization-group method may be 
regarded as a kind of variational method. In many cases, the variational calculation yields worse 
estimation for the correlation functions. 

The correlations Cpy B (#,21) and 0* tTing (9, 21) are plotted in Fig. 5 (a) and (b), respectively. 
The maxima of the correlations are located around 9 — |. This is also reported in the previous 
subsection for the system with J' —>■ oo. We hence observe that the most relevant angle for the 
correlations is given by 9 = it/ 2 as is expected from the discussions in §2. In order to estimate the 
infinite distance limit of the correlations, we plotted Cpy B (|, \i — j\) and 0* tTins (~ > N — j\) against 
l/\i — j\ in Fig. 6 (a) and (b), respectively. The correlation (||) remains finite; 

O z mB (^-,\i-j\^oo)^0M. (10) 



2' 

(The correlation would have developed fully as O z (^) = 1, if the ground state satisfied the 
condition (§) completely.) On the other hand, it is seen that the correlation (|6]) is very small. 
Considering the precision of the numerical calculators, we conclude that the correlation (||) is of 
short range. 

Finally, we show the results of the correlation (|5p with the angle 9 = 0. Note that this is 
nothing but the Neel correlation. We present a semi- logarithmic plot of the correlation against 
\i— j\; see Fig. 7. It is seen that it does decrease exponentially. The correlation length is somewhat 
consistent with an estimate in Ref. [H. 



3.3 Mechanism of the RVB order with 9 = | 

In this subsection, we clarify the essential mechanism of the development of the correlation (|J) 
with 9 = | in the region J' > 0. 

For this purpose, let us decouple the four-leg ladder system into two ladder systems. The 
Hamiltonian is given by 

4 L 

7~L = ' 

1=1 1=1 

L 

i=l 
L 

i=l 
L 

+j'J2s 3 , l -s,, l . (ii) 

i=i 

The effect of the inter-ladder coupling J" is analysed with the aid of numerical simulations below. 
Reigrotzki et al.^KQ suggested that the decoupled system (J" = 0) is a convenient starting point 
for investigating the four-leg ladder. 

First, we consider the completely decoupled case J" = 0. The hidden correlation (|5|) is reduced 

to 

Oi YB (9,\i-j\) = 2Ke ie n:^^)(e ie n:^) 

+ 2( e ie n= !+ i °* a *) ( e ie IX- ^V/), (12) 

where of = Sf i + and t* = + S^ i+ i. The brackets here denote the ground-state 

expectation value of a single two-leg ladder system. 



In Fig. 8, we show an example of spin configurations {of} and {t? } for each single-ladder 
system. The configurations are expressed in terms of the height of steps. The upward (downward) 
step at the position i of the upper drawing stands for of = 1 (of = — 1). The same rule applies 
to the lower drawing. The example here satisfies the conditions 











< 1 and 




k=i 




k=i 



< 1 



(13) 



for v z, v j 



therefore the sites with the magnetization 1 and —1 appear alternately with the sites 
with being inserted between them. Indeed, such configurations are expected to dominate the 



ground state of the two- leg ladder system (TT|, |T^, [13] , according to the RVB theory for the two- leg 
ladder. RVB pattern as in Fig. 9 satisfies the condition (13). We calculated the two terms in eq. 
flT2| ) for the exemplified configuration in Fig. 8 fixing the angle 9 to 9 = |. The expectation values 
indicated in Fig. 8 are the average over the depicted configuration and the reflected (of — > —of, 
rf — > — rf ) one. We can see that the correlation ( |T2"D does not vanish only at the sites where the 
steps of the upper and lower ladders synchronize coincidently. 

We can also show after a similar analysis that the hidden correlation (|5]) does not develop for 

9 = and 9 = tt, if we restrict the configurations with the condition ([13]). This is the same as 
what we observe in the previous subsection. Thus, the analysis based on the decoupled system 
appears to be fairly relevant. 

Now we investigate the effect of the inter-ladder coupling in terms of the above RVB picture 
numerically. As a consequence, we present more detailed information on the RVB pattern. In 
Fig. 10, we show the hidden correlation @ for the system with L = 6, J' = 1 and J" varied. We 
observe that the inter-ladder coupling J" stabilizes the hidden correlation. Recall that the essential 
mechanism of the development of the RVB correlation is the synchronization of the configurations 
of the upper-half ladder and the lower-half ladder; see Fig. 8. The RVB correlation C|y B for 
J" = expresses the contribution of the coincidental synchronization. We can observe from Fig. 

10 that the inter-ladder coupling enhances the synchronization. 

We can explain the reason of the enhancement as follows. According to the RVB picture, the 
ground state is stabilized by the resonance between the configuration of two vertical singlets and 
that of two horizontal singlets. |9|, [11]] When the inter-ladder coupling J" is turned on, the resonance 
shown in Fig. 11 becomes possible. To take advantage of the resonance, singlets along the upper- 
half ladder and those along the lower-half ladder may tend to appear at the same position i. This 
effect can enhance the synchronization of the configurations. 



4 Criticality and the phase diagram 

In this section, we investigate the development of the hidden correlations (|5|) and @ and their 
criticality when we change the parameter J'. We show the results of the exact-diagonalization 
method under the periodic-boundary condition. The result is consistent with that of the scaling 
theory [0]. 

We plotted the values of 0^y B (^, |f) and 0* tTing (~, -|) for the system of the size L = 4 and 6 
in Fig. 12. As is expected, we observe that the correlation Crvb(|) develops in the region J' > 0, 
while the correlation 0| tri (|) develops in the region J' < 0. We readily know the values of the 
correlations in two limiting cases. In the limit J' — ► oo, we showed Orvb(1) = 0.386895. In the 
limit J' — > — oo, the corelation (|J) corresponds to the string correlation of the S = 2 Heisenberg 
chain. It is reported [[24], |25[] that we have Coring (f ) ~ 

In order to see the development of the correlations and the criticality, we calculated the squared 
order parameters, 

(0U(|)0rvb(|)), (14) 



and ^ 

( ^string ( 2 ) ^string ( 2') )> (15) 

where 

Orvb(^) = 7 E e^^^+^'^Hi)^ + s» m + s* + <?I, i+1 ) (16) 
^ i=i 

and 

<W*) = j E e^Sti^ + ^ + ^ + ^(^ i + S£, + S*< + sy . (17) 
L 1=1 

We plotted the results in Fig. 13 (a) and (b), scaling the data by the factor L*. This is because 
of the following reason. At the critical point J' c = 0, the correlations are reduced to 

O z a (9,r) = (O z s=1/2 (9,r)y (a = string, RVB) (18) 



V 

'' B =x/*V>\'--J\) = \°i* ™ ^ 
If we know the correlation exponent of (|I5 ) defined in 



where C| =1 / 2 denotes the string correlation for the S 1 = 1/2 Heisenberg chain: [[37 



e>!=i /2 (0, N-il) = (s;<PZ£*is?) . (19) 



^ =1/2 (^,r)~r-^W, (20) 

the relation fllTf ) is immediately followed by 

Va (0) = 4r} S=1/2 (9) (21) 

both for a = string and a = RVB. Generalizing Hida's analysis|37] we obtain the formula 
Vs=i/2(9) = |(- ) 2 . Therefore, the exponent is given by ^ a (|) = |. 

The scaled order parameters are invariant at the critical point for various system sizes. In 
Fig. 13 (a) and (b), we observe the crossing points close to the expected critical point J' c = O.flTT) 



(Another crossing point which appears in Fig. 13 (b) around J' w 1.5 may have arose owing to 
the limited system size.) The hidden correlations appear to obey the behavior explained in §2. 



5 Summary 

The ground state of the four-chain ladder model with both ferromagnetic and antiferromagnetic 
inter-chain coupling J' has been analysed by means of the hidden correlations @ and @. The 
hidden correlations can indicate the development of the corresponding RVB states. 

The hidden correlation @ develops in the phase J' > 0, while the hidden correlation @ 
develops in the phase J' < 0. Though both the two phases are disordered, they are characterized 
in terms of the hidden long-range orders. It indicates that the RVB structure varies drastically 
at the critical point J' c = 0. These results are consistent with a recent proposal |TD[ that the point 



J' — ► oo (— oo) is the fixed point of the region J' > (< 0). A detailed analysis of RVB pattern 
has been reported for the phase J' > 0. It implies that valence bonds formed along the ladder 



are not negligible as are in the case of the double-chain ladder model. [JTTJ] This tendency may be 
more apparent as the number of the legs is increased. 

We are in position to conclude that, in the phase J' > 0, the ground state of the ladder models 
with two and four legs is the RVB state proposed by White et al. 0. The RVB pattern may be 
common to the arbitrary ladder models with even legs. 
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Figure captions 

Fig. 1 Schematic drawings of the expected RVB states. The RVB pattern of the type (a) ((b)) 
is dominant in the phase J' > (J' < 0). 

Fig. 2 Examples of spin configurations. The configurations (a) and (b) are arranged so as to 
satisfy the RVB patterns in Fig. 1 (a) and (b), respectively. The composite magnetization Sf is 
shown below each configurations. The full string order develops in both cases. 

Fig. 3 The long-range hidden correlation 0§ {VB (9, — > oo) against the angle 9 at J' — > oo. 

Fig. 4 The relative error of the ground-state energy that is calculated by means of the density- 
matrix renormalizat ion-group method for the system with the parameter J' = 1 and the system 
size L = 6 under the open-boundary condition. 

Fig. 5 The hidden correlations (a) 0py B (0, 21) and (b) O^ tTing (0, 21) with the angle 8 varied. 
The maxima are located around 9 — | . 

Fig. 6 The hidden correlations (a) ©Ryed > N — j\) an d (b) Coring (f > N ~ J\) plotted against 
the inverse of the distance — j\. 

Fig. 7 A semi-logarithmic plot of C^ VB (0, \i ~ j\) against \i — j\. For reference, we show the 
expected long-range behavior oc e-l^'l/ 5 - 5 !! as the broken line. 

Fig. 8 An example of the spin configuration for the upper double-chain ladder and the lower 
ladder. Two terms that appear in eq. (|P2|) are calculated for the configurations, respectively. 
The expectation values are the average over the depicted configuration and the reflected one 
« -> -of, 7f -> -rf). 

Fig. 9 Example of the expected RVB pattern for the decoupled (J" = 0) system. 

Fig. 10 The hidden correlation 0rvb(| > f) f° r ^ ne system with L — 6, J' — 1 and J" varied. 

Fig. 11 As the inter- ladder coupling J" is turned on, the resonance between the upper and 
lower ladders appears. 

Fig. 12 The J' dependence of the hidden correlations Crvb(| ' § ) an< ^ ^string (f ? f )■ 

Fig. 13 The J' dependance of the scaled hidden orders: (a) (Obvb(| )^RVB(f )) and (b) 

( ^string ( § ) ^string ( f ))■ The intersection point may indicate the critical point. 
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